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Abstract 

As pointed out by Lingenfelter and Ramaty (1977), the shapes of some 7-ray lines 
produced by cosmic-ray interactions with the interstellar medium potentially con- 
tain valuable information on the physical properties of dust grains, including their 
compositions and size distributions. The most promising of such lines are at 847, 
1369, 1779 and 6129 keV, from 56 Fe*, 24 Mg*, 28 Si* and 16 0*, respectively. We 
performed detailed calculations of their profiles using, in particular, available lab- 
oratory measurements combined with optical model calculations to evaluate the 
energy distributions of the recoiling excited nuclei. We show that the line shapes 
are mainly sensitive to relatively large interstellar grains, with radii ^ 0.25 fim. 
Line fluxes from the inner Galaxy are then predicted. 
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1 Introduction 



Observations of nuclear interaction 7-ray lines from the interstellar medium 
(ISM) would provide a unique tool to study Galactic cosmic-ray ions at non- 
relativistic energies, as well as the physical conditions of the emitting regions. 
If the lines produced in the gaseous phase are expected to be significantly 
Doppler-broadened, some lines produced in interstellar dust grains can be 
very narrow, because s ome of the excited nuclei can stop in solid material s 



before emitting 7-rays f Lingenfelter and Ramatvl Il977t iRamatv et a l.. 1979) 



The latter are prime candidates for detection with 7-ray telescopes having 
high spectral resolution, such as the INTEGRAL spectrometer (SPI). 

An illustrative 7-ray spectrum is shown in figure (1). It was obtained in a re- 
cent experiment performed at Orsay (Kiener et al., in preparation), in which 
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Fig. 1. Observed 7-ray spectrum from the bombardment of a thick sample of the 
Allende meteorite with 10-MeV protons. The most intense lines are labeled with 
their nominal energies and the target nuclei from which the 7-rays are produced. 
For the 6129-keV line, the single and double escape peaks are also indicated. 

a thick sample of the Allende meteorite (~5 mm in diameter) was bombarded 
with 10-MeV protons. The Allende meteorite belongs to the class of carbona- 
ceous chondrites and has a composition close to cosmic, except for the volatile 
elements. The strongest lines (apart from the 511-keV emission) originate from 
reactions with the most abundant isotopes: 16 0, 24 Mg, 28 Si and 56 Fe. One point 
of particular interest is that the 16 0* line at 6129 keV is much narrower than 
the two other 16 0* lines at 6916 and 7115 keV. This is because the 6.13-MeV 
state is relatively long-lived (Ti/ 2 =18.4 ps), such that the recoiling 16 0* nuclei 
can come to rest in the target before the 6.129-keV 7-ray is emitted, whereas 
the 6.92- and 7.12-MeV states (Ti/ 2 =4.7 and 8.3 fs, respectively) mostly de- 
excite in flight. The 6129-keV line, as well as the intense lines at 847 keV from 
56 Fe* (21/2=6.1 ps), 1369 keV from 24 Mg* (T 1/2 =1.35 ps) and 1779 keV from 
28 Si* (Ti/2=475 fs) are the prime candidates to scrutinize the interstellar dust 
grains through their 7-ray emission. 

We have calculated in detail the shapes of these four lines as they are produced 
by cosmic-ray interactions with the ISM. The interaction model and the results 
are presented in sections (2) and (3), respectively. In section (4), we evaluate 
the intensities of the predicted emissions from the inner Galaxy. 



2 Interaction model 



To calculate the line profiles, we used a Monte-Carlo method similar to that 
described in Ramatv et al. ( 1979f b Each event in the simulation corresponds 



2 




Fig. 2. Differential inelastic scattering cross sections for the 
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Mg(p,p'Y 4 Mg* l3m at E p =35.2 MeV and (b) 
Solid curves - coupled-channel calculations. 
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29.81 MeV. 



to a 7-ray producing nuclear reaction of a cosmic-ray proton or a-particle with 
a target nucleus either in the gas or locked up in a spherical and homogeneous 
dust grain of the ISM. The 7-rays can be produced by inelastic scattering reac- 
tions, e.g. 16 0(p,p'7 6 i29) 16 0, by (p, n) reactions followed by electron capture or 
positron emission, e.g. 
spallation reactio ns, e. 
iKozlovskv et al. 



56 




Fe(p,n) 56 Co(e-/5 + ) 56 Fe* 47 (T 1/2 ( 56 Co)=77.2 d), or by 
20 Ne(p,pa7 6 i 29 ) 16 O. The total cross sections are from 



except for the reactions 24 Mg(p, n) 24 Al(/? + ) 24 Mg 



1369 



=270.3 ms . We 



(T 1/2 ( 24 A1)=2.053 s) and 28 Si(p, n) 28 P(/3+) 28 Si^ 779 (T 1/2 ( 28 P) 
estimated the con tribution of the 24 Mg(p, n) 24 Al channel from the data of 
iKiang et al.l (|l989f ) and assumed the same relative contribution for the reac- 
tion 28 Si(p,n) 28 P. 

The energy distribution of the recoiling excited nuclei were calculated from 
the differential cross sections. The latter were obtained from a large number 
of experimental data, c ompleted with extensive coupled-channel calculations 
with the code ECIS94 (|Ravnall 1994). Two examples of differential inelastic 
scattering cross sections are shown in figure (2). 



The stopping powers of the nuclei r e coilin g in the grain material were calcu- 
lated with the code TRIM (|Ziederl . Il98 5). For these calculations, the grain 
composition was ass umed to be (MgSi Fe)Q/t, which is characteristic of inter- 



stellar silicates (e.g. Dwek et al. . 19971 ). 



To estimate the cosmic-r ay proton in t erstel lar spectrum, we used the disk- 
halo propagation mode l ofldones et alJ (l200ll). Coulomb a nd ionization energy 
losses were taken from iMannheim and Schlickeiserl (J1994J). For the proton in- 



elastic cross sections, we used the empirical formula given in lMoskalenko et al 



Solar-modulated spectra, calculated from the force-field approxima- 
tion, were fitted to measured proton fluxes (fig. 3). We found the best fit to 
be provided by the source spectrum Q{p)dp cx p~ 2A dp, where p is the pro- 
ton momentum. For simplicity, we used the same form for the a-particle LIS 
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Fig. 3. Cosmic-ray proton local interstellar spectrum (LIS; solid curve) and so- 
lar-modulated spectra for two forc e-field p otentials, d>=550 an d 750 MV (dashed 
curves). AMS: lAlcaraz et al.l (1200(1 ): IMAX: lMenn et al.1 fcOQOh . 



spectrum as for the protons, with an abundance ratio a/p=0.1. 



3 Gamma-ray line profiles 



Calculated 7-ray line profiles are shown in figure (4) for different interstellar 
grain size distributions. We assumed that all the available refractory elements 
Mg, Si and Fe, and half of the O are locked up in silicate grains, whereas 
the volatile elements Ne and S (which contribute through spallation reactions 
to the 6129- and 1779-keV lines, respec t ively) are in the gaseous phase (see 
Savage and Sembachl . Il99(4 iDwek et"aL ,_ 1997^1 . The do tted curves were ob- 
tained for the commonly used MRN ( Mathis et al!ll977 ) size distribution : an 



a" 3 ' 5 power law in grain radii a from a m j n =5 nm (see IDwek et all Il997l ) to 
a max =0.25 /iiii. We see that there are little differences between these spectra 
and those calculated assuming that all the target nuclei are in the interstellar 
gas (thick, grey curves). In particular, the two profiles of the 6129-keV line 
are nearly identical, because almost all of the 16 O nuclei excited within dust 
grains of radius <0.25 /im escape from the solids before decaying. 



The dashed curves were calculated for the same power law distribution, but 
with a max =2 nm. The extension of the grain radii beyond the MRN limit 
of 0.25 ixm is motivated by the observations with dust detectors aboard the 
Ulysses and Galil eo spacecrafts of relati vely large interstellar grains entering 
the solar system ( Landeraf et al. . 2000l). as well as b y models of scattering 
halos observed around X-ray sources ([Witt et all 120011 1 . The thin, solid curves 
were obtained by further extending the MRN distribution to a max =10 /im. 
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Fig. 4. Profiles of the 7-ray lines at 847 keV ( 56 Fe*), 1369 keV ( 24 Mg*), 1779 keV 
( 28 Si*) and 6129 keV ( 16 0*), excited in cosmic-ray proton and a-particle interac- 
tions with interstellar gas and dust grains. The calculated spectra are normalized 
to one photon emitted in each line and convolved with the SPI response function. 
Thin curves - grain size distributions following an a -3,5 power law in radii from 
a m j n =5 nm to a max =0.25 //m (dotted curves), 2 /xm (dashed curves) and 10 ^m 
(solid curves). For each 7-ray line, the fractions of target nuclei assumed to be in 
the grains are indicated. Thick, grey curves - all target nuclei are assumed to be in 
the gas. 

Such large grains are typ ical of presolar grains of Stardust found in primitive 
meteorites ( ZinnerL Il998l ) and should exist in various circumstellar media. We 
see that for these two grain size distributions, the very narrow component 
of the 24 Mg*, 28 Si* and 16 0* lines could in principle be resolved with SPI. 
The sensitivity of these line shapes to micrometer-sized particles provides a 
promising method for tracing large dust grains in the ISM. 



4 Gamma-ray line fluxes from the inner Galaxy 



The 7-ray line intensities from the inner Galaxy were estimated by nor- 
malizing the emissivity calculations to the observed flux of high-energy 7- 
rays (>70 MeV) du e to 7r ° production and decay. We derived the latter 
from iBiisching et al.l ([200lL fig. 1; EGRET data of observing phases 1-4): 
$^°~8.3xl0- 5 photons cm" 2 s" 1 , for -40°<£<40° and -6°<6<6°. We used 
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the model of Dermei (1986) to calculate the production of 7r°-decay 7-rays by 
cosmic-rays with the LIS spectrum of figure (3). We assumed that the metal 
abundances in the inner Galaxy are on the average twice solar. We then ob- 
tained for the fluxes of the 847-, 1369-, 1779- and 6129-keV lines excited in 
proton and a-particle interactions with both interstellar gas and dust grains: 
5.6xlCT 8 , 7.5 xlCT 8 , 4.3 xlCT 8 and 3.1 xlCT 7 photons cm -2 s -1 , respectively. 
In comparison, the calculated flux of the relatively strong line at 4438 keV 
from 12 C* (FWHM~150 keV) is 7.1xl(T 7 photons cm" 2 s" 1 and that of the 
so-called a— a line at ~450 keV (complex line shape with FWHM~80 keV; 
see Tatischeff et al., 2001) is 3.9xl0~ 7 photons cm -2 s _1 . All these diffuse 
emission fluxes are far below the SPI sensitivity, such that a near future de- 
tection of nuclear interaction 7-ray lines from the inner Galaxy is unlikely, 
unless there is a large population of Galactic cosmic-rays with kinetic energies 
below the threshold for tt° production (290 MeV for p+H collisions). Such a 
distinct cosmic-ray component predominant at low energies has been proposed 
to account for the quasi line ar correlation betw een Be and Fe abundances for 
metallicities [Fe/H]<-1 (e.g. Casse et al. . 19951 ). 
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